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We study a model of two interacting levels that are attached to two electronic leads, where one 
of the levels is attached very weakly to the leads. We use rate equations to calculate the average 
current and the noise of electrons transmitted through the two levels. We show that the shot noise 
is enhanced because of the interactions and that the Fano factor depends on the properties of the 
couplings between the levels and the leads. We study both sequential tunneling and cotunneling 
processes and show that there is a range of parameters in which the cotunneling processes affect the 
noise significantly, even though most of the current is carried by sequential tunneling processes. 
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I. INTRODUCTION 

Fluctuations in electrical current, which we simply call 
noise provide additional information about the charge 
transport that is not accessible from average current mea- 
surements (for a review see Ref. Among the various 
noise sources we focus on the shot noise. The discrete- 
ness of the transferred charge causes fluctuations in the 
current. These fluctuations, named shot noise, depend 
on the charge of the conducting particles and therefore 
measurements of shot noise provide information on the 
discrete nature of the conducting particles and their cor- 
relations. The shot noise of a Poisson process of un- 
correlated current pulses of charge e is ^shot = 2e(/), 
where (/) is the time averaged current. The noise is 
proportional to the average current since in a Poisson 
distribution the mean equals the variance, therefore the 
mean number of current pulses equals its variance. The 
Fano factor, F = §f^j, is a dimensionless parameter that 
characterizes the granularity of the current. When driv- 
ing a current through a single spinless electronic level, 
the correlations that are imposed by the Pauli exclusion 
principle reduce the Fano factor; the Fano factor varies 
between half and one depending on the symmetry of the 
couplings between the level and the external leads^Ti^. 

Correlations can also be imposed by Coulomb 
interactions^. In most cases, the effect of Coulomb re- 
pulsion on the noise of a mesoscopic system is similar 
to the effect of the Pauli exclusion principle. Both im- 
pose a time delay between consecutive current pulses 
and therefore we expect negative correlations between 
them that suppress the shot noised. However, several 
theoretical works on various systems have showed that 
Coulomb interactions might also lead to a super-Poisson 
shot noise with a Fano factor F > 1. Examples of such 
systems are quantum dots that are coupled to ferromag- 
netic leads^r— , multi-levels quantum dot o^^'^^ , multi- 
dots structures'^"—, and also three terminal quantum 
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There are also experimental works 



which a super-Poisson noise was measured in quantum 
dots, rather than the sub-Poisson noise, which is expected 
from the single level model of the quantum dots. 



A simple mechanism that might explain the enhance- 
ment of the Fano factor in such systems is tunneling 
through two levels that are coupled to the leads, where 
one of the levels is coupled much stronger than the other 
level. The two levels are interacting, namely, there is a 
Coulomb repulsion between electrons that occupy the two 
levels. In this case, the electrons that cross the system 
tunnel mainly through the level that is strongly coupled 
to the leads. However, once in a while an electron can 
tunnel into the weakly coupled level and then, because 
of the Coulomb interactions, the tunneling of other elec- 
trons through the strongly coupled level is prevented, and 
the current is blocked. The current resumes only after the 
electron tunnels out of the weakly coupled level. There- 
fore the intuitive picture is a current that is blocked oc- 
casionally and therefore the noise is enhanced. The idea 
of two interacting levels as a possible source of super- 
Poisson noise was discussed in the context of quantum 
dots that are coupled to ferromagnetic leads^i^ and also 
in the context of double quantum dots structures^. 

In this work, we analyze in detail the model of two lev- 
els with Coulomb interactions that are attached to two 
electronic leads. In particular, we study how the shot 
noise and the Fano factor depend on the left-right asym- 
metry, namely, the asymmetry between the couplings to 
the two external leads. The fact that the enhancement of 
the shot noise depends on the left-right asymmetry of the 
coupling to the leads emerges from previous works (e.g., 
the results of Refs. [l^ and fisl). nevertheless a complete 
theoretical analysis of this dependence is missing. We 
find that the Coulomb interactions affect significantly the 
dependence of the Fano factor on the left-right asymme- 
try. We study how the temperature affects this depen- 
dence and show that finite temperature suppresses the 
Fano factor at very asymmetric couplings. In addition, 
we include in our calculations not only the sequential 
tunneling processes (leading order perturbation in the 
tunneling coefficients) but also the cotunneling processes 
(next leading order perturbation) and we find that the 
latter play an important role in this system. Surprisingly, 
and perhaps this is the main point of our work, there is 
a range of parameters where the noise is governed by the 
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rare cotunneling processes even though most of the aver- 
age current is carried by sequential tunnehng processes. 
The importance of the cotunnehng processes follows the 
fact that they allow the system to change its occupation 
from one level to the other level in a single quantum pro- 
cess. The two levels mechanism for noise enhancement 
can be found in many physical realizations such as single 
level with spin dependant coupling and double quantum 
dots or two levels in a quantum dot in a strong magnetic 
field. Although the quantitative details of each system 
are different, the qualitative behavior of the noise en- 
hancement is the same. 

The paper is organized as follows: in Sec. |TT1 we present 
the model that we analyze, then we give in Sec. Illll a sim- 
plified qualitative description of the transport through 
the two levels, which we use to develop an intuitive un- 
derstanding of the results that we obtain later through 
the rate equations method. Next, in Sec. IIVI we give 
the details of the theoretical calculations. In Sec. |Vl we 
present the results of our calculations and discuss them. 
Last, we summarize our results in Sec. IVII In the two ap- 
pendices, we list all the tunneling rates that are relevant 
for the theoretical calculations (Appendix |X| and give 
technical details about the regularization procedure that 
we use in the calculations of the cotunneling processes' 
rates (Appendix |B|) . 



II. MODEL 

We consider a model of two interacting levels that are 
attached to two leads (see Fig. [T]). For simplicity, we 
discuss a spinless problem^l. One of the levels (level 2 in 
Fig. [T]) is attached very weakly (compared to the other 
level) to the leads. A simultaneous occupation of the 
two levels is possible, however this situation is not likely 
to happen as it requires an additional charging energy, 
U, because of the Coulomb interaction. This model de- 
scribes, for example, an interacting two level quantum 
dot in a strong magnetic field. The Hamiltonian that 
describes the system is 



H = H\ 



leads 
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(1) 



where 



Eid\di + E2dld2 + Ud\did\d2, 

H.c. 



^,t^Lld, + tfRld 

i.k 



Here, Lk (Rk) are left (right) lead annihilation opera- 
tors, di is the ith level annihilation operator, and we 
have assumed that the tunneling coefficients {tf'^) are 
independent of the energy. We also assume that all the 
tunneling coefficients have the same sign. 

We calculate the average current and the zero- 
frequency noise using rate equations method in which 
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Pi=7r,.{\tf\^ + \tf\') 
^l2 = ^H\t^\' + \t2\') 



Figure 1: A schematic picture of the model. The energies of 
the two levels are Ei and E2 . The bias voltage is the difference 
between the two chemical potentials: eVbias = fJ'^ — M^- Here 
the width of level 2 is much smaller than the width of level 1 
which is assumed to be much smaller than the bias voltages: 
72 ^ 71 ^ IfJ-^'^ — £'1,2 1 . We assume that the tunneling 
coefficients are independent of the energy and we also assume 
for simplicity that they all have the same sign, and that the 
two ratios between the left tunneling coefficients and the right 



tunneling coefficients are identical, ti /t^ = t2 /t 



we assume that the two levels are weakly coupled to the 
leads. By weakly coupled we mean that the widths of 
the levels are much smaller than the temperature, or al- 
ternatively, at low temperatures it means that the bias 
voltage (see Fig. [T]) is much larger than the levels' widths 



R.L 



E, 



(2) 



where /i^'^ are the electro-chemical potentials of the 
leads and we have assumed for simplicity the same den- 
sity of states, V, in the two leads. 



III. QUALITATIVE SIMPLIFIED MODEL 

In this section, we study a simplified intuitive model 
that captures, at least qualitatively, most of the results 
that we later achieve through a more rigorous analysis. 
The propose of this section is to establish a simple phys- 
ical picture that we can use to interpret the results that 
we get through the rate equations formalism. Consider 
the two level system that was discussed in the previous 
section and depicted in Fig. [T] and assume strong inter- 
actions, i.e., a large U. At this point, we also assume, 
for simplicity, zero temperature. As a function of Ei, 
there are two regions where we expect to have current: 
at < El < fJ" where the two levels are between the 
chemical potentials (assuming that \Ei — i?2| ^ Vbias), 
and at < Ei + U < fi^ where one of the levels is occu- 
pied and effectively, because of the Coulomb interaction, 
the other level is shifted up and placed between the chem- 
ical potentials. We choose to focus in this section on the 
later: fi^ <Ei + U < fi^. 
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Since level 2 is coupled weaker than level 1 to the leads, 
most of the time the current flows through level 1 {i.e. 
level 2 is occupied and electrons enter and leave level 
1). However, after a while, the electron in level 2 can 
tunnel out to the right lead and by that, because of the 
Coulomb interaction, it reduces the effective energy of 
level 1 (from Ei + U to Ei) making the tunneling out 
of level 1 impossible. The current is therefore blocked. 
The current resumes only after a new electron from the 
left lead tunnels into level 2. Hence, the picture is a 
current (through level 1) that is stopped occasionally (by 
tunneling out of level 2). This situation is schematically 
drawn in Fig. [2] 

A. Simplified model: shot noise on top of a 
telegraph noise 

The current, as drawn in Fig. [21 fluctuates between 
two modes: A zero-current mode where the current is 
dramatically suppressed, and a nonzero mode where the 
current is carried by pulses of charge (electrons) that tun- 
nel through level 1. We therefore suggest the following 
simplified model: the current I(t) is a multiplication of 
two signals 

graph (3) 

where /shot(0 the current through level 1 and 
Ctciograph(i) is a random telegraph signal that fluctuates 
between two values, zero and one, according to the oc- 
cupation of level 2. In this simplified model, the current 
flows only through level 1 and only when level 2 is oc- 
cupied. We neglect the effects that the Pauli principle 
imposes on the current /shot (t) and treat it as a sequence 
of current pulses of charge e with a Poissonian statistics 
characterized by a rate Fi. The changes in the value of 
Ctcicgraph are Poissonian events with the rates 1/ti and 
I/tq for the events of changing the value of Ctcicgraph 
from one to zero and zero to one, respectively. The fluc- 
tuations in the current /shot(i) are known as shot noise, 
and the fluctuations in the signal Ctcicgraph (i) are known 
as telegraph noisoS^. 

Assuming that /shot and Ctcicgraph are uncorrelated, 
the average current is 

(/) = (/shot) (Ctcicgraph) = eFi ^ . (4) 

Tl + To 

The noise is related to the auto-correlation function of 
the total current via Wiener-Khinchin theorem 

/oo 
rfre-((/(i + r)/(i)>-(/>^) (5) 
-oo 

= (/shot)^'5'tclcgraph('^) + (Ctcicgraph) ^ "^shot (w) 

+ — S'shot(l^) * 'S'tclcgraph(w)- 
TT 

Using the known results for the shot noise and the tele- 
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Level '2' is occupied 




Level '2' is empty 

Figure 2: A qualitative description of the transport through 
the two level system in the region jj.^ < E\ + U < Elec- 
trons are tunneling through level 1 until the electron in level 
2 tunnels out and then the current is blocked. The tunneling 
events through level 1 resume when another electron tunnels 
into level 2. The current behaves as a multiplication of two 
signals; A sequence of current pulses through level 1 and a 
telegraphic signal, the occupation of level 2. Fi is the aver- 
age rate of tunneling events through level 1, ro is the average 
time duration until level 2 is filled, and ri is the average time 
duration until level 2 gets empty. 

graph noiseiSSi we get 

^'^^ ^ Vl+To (Tl+To)((l/Tl + l/ro)2+w2)• 

(6) 



B. Fano factor and asymmetry dependence 

Before we study the Fano factor and its left-right asym- 
metry dependence, it seems necessary to inquire about 
the linear dependence of the noise on the current as the 
second term in ^ is quadratic in Fi. Since the rates 
1/ti and I/tq are the rates of tunneling out of and into 
level 2, they depend, up to a symmetry factor, linearly 
on 72 that was defined in Eq. Similarly, Fi is linear 
in 7i. Therefore, under the reasonable assumption that 
71/72 is independent of the left-right asymmetry (that is, 
the couplings of the levels to the leads cannot be changed 
independently), 1/ti and 1/to are linear in Fi. Hence, 
the zero- frequency noise, 5(0), depends linearly on Fi 
and therefore depends linearly on the current. 

The probability of finding level 2 occupied, which is 
the probability of finding Ctcicgraph = 1, is Pf = 
Similarly, the probability of finding level 2 empty, which 
at zero temperature is the probability of finding only level 
1 occupied is = — . The Fano factor is 

Notice that FiTi is the average number of tunneling 
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events during a Ctciograph = 1 stage, which is the av- 
erage number of tunneUng events through level 1 before 
a tunneling event out of level 2 takes place. At zero 
temperature it can be estimated as |tf P/I^I^P, or, as- 
suming the same left-right asymmetry for the two levels. 
Tin =71/72- The Fano factor ([7]) becomes 

F = 1 + 2— (8) 

72 

The probability depends on the asymmetry between 
the coupling to the left and the right leads, i.e. on the 
ratio |if p. To see this consider the simple case of 
zero temperature and large U. For /i^ > /i^ (see Fig. 
[T|), the rate of tunneling out of level 2 is, according to 
Fermi's golden rule, Ff^o = x 1*2^1^- Similarly, the rate 
of tunneling into level 2 is Fe->.f = xl*2 If we neglect 
cotunneling effects there are no direct tunneling processes 
from level 1 to level 2, so the steady state probabilities 
satisfy 

Pcrc^f = PfFf^e. (9) 

Pf is the probability of finding both level 1 and level 2 
occupied. In the limit \Ei — i?2| ^ Vbias, the probability 
of finding only level 1 occupied, Pe, equals the probability 
of finding only level 2 occupied. In the limit of large U, 
the probability of finding both level 1 and level 2 empty 
(for < El + U < fi^) is zero. Thus, Pf + 2Pc = 1. 
The steady state solution ([9]) becomes 

By increasing |tf we decrease Pq. In the limit 

l^f lVl*f P ^ 00, Po ^ so the Fano factor © P 1. 
By decreasing \tf\'^/\tf'\'^ we increase Pc and in the limit 
|i,f P/|tf P it gets its maximal value Pc = 1/2 and 

the Fano factor is maximal. P = 1 + ^ w t?^. 

272 272 



and forth between the right lead and the system before 
a tunneling event from the left lead to the system takes 
place. The time scales ti and tq in this case, become 
smaller than l/Fi, and the value of FiTi goes to zero. 
Therefore, we expect that at finite temperature the Fano 
factor will have the value one in the limit |if P/|tf P ^ 0. 
To conclude, as we decrease |tf P/|tf P the Fano factor 
([7|) gets larger, but at finite temperature, at some point, if 
we decrease l^fP/lif P even further, the Fano factor will 
start to decrease toward the value one at |if P/|tf P — > 0. 



2. Cotunneling effects 

Similar to finite temperature, cotunneling processes 
may suppress the Fano factor at |if P/|tf P — 0. Con- 
sider the limit |tf p < |tf p. Without cotunnehng, if 
level 2 is empty, the electron in level 1 need to wait a 
long time before it can tunnel out to the right lead since 
such a tunneling event must follows a tunneling event 
into level 2 which, at zero temperature, is possible only 
from the left lead. Cotunneling processes, however, allow 
the two processes at once; occupying level 2 and evacu- 
ating level 1 in a single quantum process. In particular, 
the electron in level 1 can virtually tunnel out to the 
right lead while another electron is virtually tunnel from 
the right lead into level 2 (we use the term virtually to 
emphasize the fact that the intermediate state does not 
conserve energy). In the limit |tf P/|if P the to- 
tal rate of such processes may become larger than the 
rate of sequential tunneling from the left lead into level 
1 . Thus the occupation of level 2 and therefore the tele- 
graphic signal, Ctciograph (i), fluctuates much faster than 
the pulses' rate 1/Fi, and the value of FiTi goes to zero. 
Therefore we expect a suppression of the Fano factor due 
to the cotunneling processes toward the value one in the 
limit ItflVltfp^O. 



1. Finite temperature 

At zero temperature, the only possible tunneling event 
that follows a tunneling event from the system to the 
right lead, is from the left lead into the empty level. 
Therefore, changing the value of Cteiegraph from one to 
zero in the simplified model, corresponds to a tunnel- 
ing event from level 2 to the right. In addition, every 
tunneling event from level 1 to the right, corresponds 
to a current-pulse in the simplified model. Hence, we 
can estimate the number of tunneling events during a 
Ctciograph = 1 stage as Fin ~ l*f P/l^fP- At finite tem- 
perature, there is also a finite probability that a tunneling 
event from the system to the right will be followed by a 
tunneling event from the right lead back to the empty 
level. In the hmit |if | Vl^f P ^ the probability of tun- 
neling from right to left may become important. If the 
couplings to the right are much stronger than the cou- 
plings to the left, electrons may tunnel many times back 



C. Weak interactions 

Up to this point, we have assumed strong interactions, 
namely, a very large U. If J7 is not large compared to 
the bias voltage, U < eVbias, changing the occupation 
of level 2 doesn't block completely the current through 
level 1 since electrons can tunnel through level 1 in both 
cases; while level 2 is empty or occupied. Yet, we can still 
use the intuitive picture of a sequence of current pulses 
through level 1 and a random telegraph signal describing 
the occupation of level 2. The rate of the pulses depends 
on the occupation of level 2 and we consider two different 
rates: Fi describes the tunneling rate through level 1 
while level 2 is full and Fi describes the tunneling rate 
through level 1 while level 2 is empty (previously Fi was 
zero ). Similar to Eq. we consider the current 

I{t) ~ -fshot(i) X Ctciograph (i)+^hot(^) X (1 ^ Ctclograph(i)) , 

(11) 
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where /shot (/shot) is a sequence of current pulses with a 
characteristic rate Fi (Fi). The Fano factor is 



F = 1 



2(Fi-fi)SVf 
(ro + ri)2(TiFi+Tofi)^ 



(12) 



where ti is the average time duration in which level 2 
is occupied and electrons tunnel through level 1 with an 
average rate Fi, and tq is the average time duration in 
which level 2 is empty and electrons tunnel through level 
1 with an average rate Fi. 

At low temperatures, KbT <C eVbias, if fJ.^ < Ei < fi^ 
and < Ei + U < jj.^ , the tunneling rate through level 1 
barely depends on the occupation of level 2, i.e. Fi « Fi. 
We can approximate by 



1 



2(AFi)- 
F1F2 



(13) 



where AFi = f 1 — Fi and F2 = (tq + ti)"-'^ (at low tem- 
peratures this is the tunneling rate through level 2). No- 
tice that Pc = ro+Ti ' "^^^icii is the probability of finding 
level 2 empty, is different from the probability of finding 
only level 1 occupied, since now the two levels can be 
empty simultaneously. 

The results of the simplified model which we presented 
in this section, can be used to analyze transport through 
other systems as well. Systems with discrete charge car- 
riers and in addition a telegraphic behavior can be an- 
alyzed in the same way. Example of such a system is a 
quantum point contact capacitively coupled to a system 
that alternates between two charge states^S. 



IV. RATE EQUATIONS FORMALISM 

The two level system is weakly coupled to the leads, 
therefore we may use the rate equations formalismi^S to 
calculate the average current and the current noise. In 
the weak coupling regime, we describe the whole sys- 
tem using the four eigenstates of the two level system: 
\ni,n2) {rii labels the eigenvalues of djdi) and treat the 
tunneling Hamiltonian, Ht, as a perturbation. Transition 
rates between these states are related to rates of electron 
tunneling events. In the following two subsections we 
first review the method of current and noise calculations 
using rate equations and we then discuss the tunneling 
processes that are relevant to our model. 



A. Rate equations 

We shortly review the technique for calculating cur- 
rent and shot noise that was developed in Ref. [30| and 
the generalization of it to include cotunneling processes 
that was developed in Ref. [3l|. Consider the two level 
system that we discuss. The system can be in one of four 
states, |ni,ri2), and its dynamics is driven by tunneling 



processes. Generally, there are a few tunneling processes 
that change the system's state from \i) to |/) which we 
denote their rates by lu"^ where a labels the different 
possible tunneling processes. The rate LUij — o;"^ is 
the total transition rate from the state \i) to the state |/). 
The system's dynamics is described by the rate equations 

^F(/, t/z) = J2 (^(^' ^A)^*^-./ - Vfc)w/.fc) , (14) 
fees 

where P{f,t/i) is the probability of the system to be in 
the state |/) at time t if it was in the state |i) at t = 0, 
so the initial condition is P{f,t = 0/i) = 5ij. Eq. 
neglects coherence superpositions of different states (the 
terms P{f,t/i) are the diagonal matrix elements of the 
system's density matrix). Neglecting the coherent super- 
positions of the states |1,0) and |0,1) is justified when 
the coherence time is much shorter than the delay time 
between consecutive tunneling events. Alternatively, if 
there is a quantum number that distinguishes the two 
states {e.g., spin), if at some point in time the density 
matrix that describes the system is diagonal, then coher- 
ent superpositions are zero at all later times. 

It appears to be useful to write the rate equations in 
a matrix form. We define the vector (there are four such 
vectors) 

p,(i) = iPioo,t/i),p{io,t/i),p{oi,t/i),p{n,t/i)), 

(15) 

and the matrix 

/ — i^nn i. Win. on uiCiiAm aJii.nn \ 



M : 



(16) 

The rate equations become ^Pj(t) = MPi{t), with 
the initial condition P,(t = 0) = (0, . . . , 0, 1, 0, . . . , 0) = 
The solution of Eq. (|14p is readily found to be 
Pi{t) — e^^^Ci. Let P** be the stationary solution, 
namely, MP^* = and (since it is a probabilities vec- 
tor)^P,f = 1. 

n 

current. We define the quantity s°^, the total number 
of electrons that tunnel through the junction a = L, R to 
the right during the process a that changes the system's 
state from |i) to |/). A negative sign corresponds to elec- 
trons that are moving to the left. In our model, sfj can 
get the values ±2, ±1,0. The stationary current through 
the junction a can be written as 
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1^01,00 


1^11,00 


k^OO 








tJOO.lO 


1^01,10 


<^11,10 




fc#10 






ti^OO.Ol 


1^10,01 
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<^00,11 


1^10,11 


1^01,11 











iJeS a 



(17) 



Zero frequency noise. The noise is related to the auto- 
correlation function via Wiener-Khinchin theorem 



Sab{uj) 



dTe'^^{{Iait + T)h{t))-{lf), (18) 



6 



where we are interested in the zero frequency hmit. We 
write here a compact expression for the noise, details 
of the derivation can be found in Ref. [3l|. The zero 
frequency noise can be written as 

SabiuJ ^ 0) = 2e2{trU,fc-Wb-M~iY,-W,-M-iY4, 

(19) 

with the foUowing vectors: 

{Vab)^ ^EE</</^'*</' (20) 

fes a 

(Ya), -EE<"^*<"^^*' (21) 

ieS a 

(Wfc)fc ^Y.Y.'''^f<f (22) 
fes a 

We have used the trace of a vector to denote the sum of 
its elements. Although the matrix M is not invertible, 
there is only one traceless vector Va {i.e. the sum of all 
its elements is zero) that satisfies MV^ = and we 

use this vector as M~^Ya. 

B. Tunneling rates 

We describe the system using four states, |rti, 712). The 
system dynamics is driven by transitions between states 
that are caused by tunneling processes. For example: 
Transition from the state |0, 0) to the state |1, 0) happens 
when an electron is tunneling from the left lead or the 
right lead into level 1 while level 2 is empty. The term 
Ht in the Hamiltonian ^ allows tunneling processes and 
the rates of the transitions are derived perturbatively in 

Sequential tunneling rates. To the lowest order in 
Ht the transition rates can be calculated using Fermi's 
golden rule. We use the notation uj~^j for the rate of a 
tunneling process that changes the system's state from 'i' 
to 'j' by tunneling an electron from the left to the right 
direction (and similarly lo^j for electron that moves from 
right to left). For example, 

^0^,10 -rf'^Fi3(i;i-A^^) (23) 

is the rate of tunneling from the left lead into level 1 
while level 2 is empty, where we have defined 

rf^.||tf^ ^f^.^|^f|^ (24) 

and the Fermi's function Fp£){x) = (1 + e^^)~^ gives 
the probability for the availability of an electron for the 
tunneling process. In Fig. |3] we depict the lowest order 
tunneling processes, to which we refer as sequential tun- 
neling processes. The rates of all the sequential tunneling 
processes are listed in Appendix \X\ 

Cotunneling rates. The next leading order perturba- 
tion in Ht generates cotunneling processes with interme- 
diate virtual states. Usually, when the tunneling coeffi- 
cients are small (71^2 of Eq. are much smaller than 
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Figure 3: A schematic picture of the lowest order in Ht transi- 
tions (the sequential tunneling processes) and their notations. 



the temperature or bias voltage) the higher order pertur- 
bation theory is not crucial as it barely improves the ap- 
proximation. However, this is not always true. Consider, 
for example, the limit ^ tf ; in this cotunnel- 
ing process, which takes an electron out of level 2 into 
the left lead and takes another electron from the left lead 
into level 1 , can become more likely to happen than a se- 
quential tunneling from level 2 to the right lead. In other 
words, a second leading order perturbation in t\ can be 
more important that a first leading order perturbation in 
t,f . We discuss later the importance of the cotunneling 
processes in the two level model and at the moment we 
emphasize that this is more than a small improvement 
of the approximation. We discuss two types of cotun- 
neling processes: elastic-cotunneling, namely, processes 
that contribute current but don't change the state of the 
two level system, and inelastic-cotunneling, namely, pro- 
cesses that change the state of the two level system (and 
might not contribute to the current through it). 

Elastic-cotunneling rates. Figure 5] depicts schemati- 
cally the elastic-cotunneling processes that we take into 
account in the transport calculations of the two level sys- 
tem. Each process has two possible intermediate states. 
For example, electron can tunnel through an empty sys- 
tem via level 1 or 2, thus, cotunneling processes of the 
form |0,0) — > |0, 0) have two possible intermediate 
states: |1,0) and |0, 1). The two possible paths interfere 
and we need to sum the amplitudes of the two possibili- 
ties rather than their probabilities. If there is a quantum 
number that distinguish the two levels ( e.g., if the system 
is a single spinful level with spin-dependent couplings to 
the leads) the two path do not interfere, and we simply 
sum their rates. We use the notation uj~l (w,^) for elastic- 
cotunneling processes in which an electron is tunneling to 
the right (left) direction (see Fig. 2]). Elastic-cotunneling 
processes in which an electron tunnels back and forth be- 
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Figure 4: A schematic picture of the elastic cotunnehng pro- 
cesses and their notations. Each process has two possible 
virtual intermediate states. The rate is the sum of the ampli- 
tudes of the two possible paths (rather than the sum of their 
probabilities). The processes cj*^ are the same as after 
exchanging the final states with the initial states. 



tween one of the levels and one of the leads don't change 
the state of the system and don't contribute any current, 
and therefore don't appear directly in the transport cal- 
culations. The total elastic-cotunneling rates are the sum 
of the rates of all the possible processes, namely, integra- 
tion over all incoming electron's energies. For example, 



''00,00 



j deFpD {e - Ai^) (l - Fpo (e - /i"")) 

(25) 



El 



The rates of all the elastic-cotunneling processes are 
listed in Appendix [Xj Equation (j25p is a formal expres- 
sion and the actual rate, which we use in the rate equa- 
tions, cannot be directly calculated from it. The rea- 
son is the divergence of the integral due to the finite 
widths of the two levels (which we treat as two delta 
functions in energy). This divergence was already dis- 
cussed before^i^, and a regularization scheme for the 
calculation of the cotunneling rates was developed. We 
summarize the regularization scheme in Appendix |B] 

It is worth noting that one can avoid the necessity of 
regularization by using the diagrammatic technique that 
was developed in Refs. (33436| . We find that additional 
correction due to levels shifts and broadening captured 
by this approach^ are irrelevant in the large bias limit 
(7i ^ eVhias) that we consider. Calculation procedure of 
the average current and current noise using this technique 
was developed in Ref. {ssj and gives the same results in 
the 7i < eFbias limit. 

Inelastic-cotunneling rates. In Fig. [51 we depict the 
cotunneling processes that change the state of the two 
level system, i.e. the inelastic-cotunneling processes. 
To this order in Ht, the inelastic-cotunneling processes 
change the system's state between the following states: 
|1,0) i — > |0, 1), |0,0) i — > |1,1). The latter are some- 
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Figure 5: A schematic picture of the inelastic cotunneling 
processes and their notations. The processes cjoi.io are the 
same as a;io,oi after exchanging the two levels. The processes 
<^ii,oo are the same as tJoo.ii after exchanging the final states 
with the initial states. 



what more complex than the other cotunneling processes 
as they have four possible intermediate states (see Fig. 
[5]). We use the notation uj^- for the rate of processes in 
which the system changes its state from \i) to |j) in the 
following way: the electron that enters or leaves level 1 
tunnels to the right direction, while the electron that en- 
ters or leaves level 2 tunnels to the left, and similarly we 
define , w^- , . We use the notation Wg^ ^ii <^n,oo 
for processes in which the two electrons enter or leave the 
two levels by tunneling one to the right and the other to 
the left. The rates of all the inelastic-cotunneling pro- 
cesses are listed in Appendix [Xj 

Beyond the cotunneling approximation. The rate equa- 
tions based calculation is valid as long as the tunneling 
coefficients are small enough as we insert the tunneling 
processes only up to second order in the perturbation 
Ht. Practically, it means that either 7i/eVbias or ^i/KsT 
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need to be small numbers. Yet, the next leading order 
in Hf generates logarithmic contributions that diverge at 
low temperatures and bias voltages^. Hence, for bias 
voltages smaller than a characteristic energy scale, the 
Kondo temperature, the perturbative approach fails. The 
Kondo temperature in our case is Tk ^ y/jiUe^^^^''^ 
and in all cases in this work we consider much larger bias 
voltages. We also want to note that small corrections due 
to the renormalization of the energy levels and broaden- 
ing play very minor role in the large bias case that we 
consider. They at most slightly modify the quantitative 
results with no important effect on the qualitative behav- 
ior. 



V. RESULTS 

In this section we present the main results for the shot 
noise and the current through the two level system that 
is modeled by Eq. ([T]). The calculations are based on 
the rate equations method which we presented in Sec. 
IIVI The results are explained using the approximated 
intuitive approach of signal analysis that we developed 

in Sec. nni 

A. Strong interactions 

The current: In Fig. Elja). we plot the symmetric cur- 
rent (rf = ) through the two level system as a func- 
tion of El at large U {U > eVbias) with the following 
parameters (F ee Tf + + + F^): KbT = hV, 
eVbias = lOhV, Ei^E2^ hV, U = lOOhT and Ff = 
2Qpfl,L Two Coulomb peaks with a width ~ eVbias ap- 
pear in the current: at i^i sa (the right peak), where 
the two levels are between the leads' chemical potentials, 
/i^ and (similar to the schematic picture in Fig. [1]), 
and at Ei « —U, where one of the levels is occupied so 
the other level is effectively shifted and placed between 
/i^ and 1^1^. The right peak is a bit lower than the left 
peak because of the small energy difference between the 
two levels {Ei > E2). At finite temperature, this small 
energy difference makes the probability of finding level 
2 occupied a bit larger than the probability of finding 
level 1 occupied. While near the right peak the current 
through the strongly coupled level- level 1- is blocked by 
the occupation of level 2, near the left peak the occupa- 
tion of level 2 allows it, and therefore the left peak is a 
bit higher. If the sign of the energy difference was the 
opposite {El < E2) the right peak was higher than the 
left peak. 

In Figs. and |6ji we plot the left-right asymmetry 
dependence of the current at the right and left peaks re- 
spectively. We assume the same left-right asymmetry for 
the two levels, Ff/Ff = F^/F^ = F^/F^. Due to the 
lack of a particle-hole symmetry at the peaks the current 
is not maximal when the system is symmetrically coupled 
to the leads. To understand this consider the current at 
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Figure 6: The current through the two level system with the 
following parameters (F = Ff -I- Ff + + F^): KbT = hV, 
eVbiaB = lOfiF, Ei-E2 = KT, U = lOOftF and Pf'-^ = lOF^'^. 
We assume the same left-right asymmetry for the two levels: 
rf /rf = r^/r^ = r^/r^. (a) The symmetric (F^/F^ = 1) 
current as a function of Ei. (b) The current at the right 
peak. El = 0, as a function of the left-right asymmetry. The 
dashed line is the sequential tunneling current, and the solid 
line is the current including both sequential and cotunnel- 
ing processes, (c) The (cotunneling) current at the valley. 
El — —U/2, as a function of the left-right asymmetry. The 
sequential tunneling current is practically zero at the valley, 
(d) The current at the left peak, Ei = —U, as a function of 
the left-right asymmetry. The dashed line is the sequential 
tunneling current, and the solid line is the current including 
both sequential and cotunneling processes. 



the left peak in the simple case of zero temperature and 
no cotunneling processes. The right lead current (fT7|) in 
this case is simply (/) = Pii(Ff -|- F^) « PnVf where 
Pii is the probability of finding the system doubly occu- 
pied. The probability Pn can be easily calculated since 
the probability of finding the system empty in this case 
(zero temperature and Ei = —U) is zero. Also, in the 
limit \Ei — E2\ ^ Vbias the probability of finding only 
level 1 occupied, Pio, and the probability of finding only 
level 2 occupied Pqi are identical and given by Pc of Eq. 
(fTU)) . The probability of finding the system doubly occu- 
pied is therefore Pu — 1 — 2Pc- The average current is 
W ~ ^1 ^ 2Pc) = (ri + -Ti ) pL /r«+i (-^ ~ r^- /r«+2 
and it is maximal, for a fixed Ff + Ff , at F^/F^ = \/2. 
Indeed, the current in Fig. [6jl is maximal at F^/F^ « -\/2 
(the calculation is done at finite temperature and in- 
cludes cotunneling therefore \/2 is only an approxima- 
tion). Similarly, the maximum of the right peak is at 
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r^/r^ « 1/a/2. The current at the valley Ei = -U/2, 
depicted in Fig. [6t, carried by cotunneling processes, is 
maximal where the two levels are symmetrically coupled 
to the leads, = T^. 

The Fano factor. The current is carried by tunneling 
of electrons, namely current pulses of charge 'e', result- 
ing in a shot noise. On top of the tunneling events, as we 
discussed in Sec. IIIIl there is also a telegraph noise; by 
tunneling into or out of level 2, we change the tunneling 
rate through level 1. Since most of the current is car- 
ried by tunneling through level 1 the current alternates 
between two different average values. For example, if we 
focus on the left peak of the current, Ei = —U (see Fig. 
EK), most of the tunneling events are via level 1 while 
level 2 is occupied. However, because of the strong inter- 
actions, each time the electron leaves the narrow level, 
level 2, the current drops dramatically and resumes only 
when a new electron enters level 2. 

n Fig. [71 we depict the Fano factor at the left peak 
{El = —U) of the current through the two level system 
with the following parameters: eVbjas = O-IU, Ei — E2 = 
O.OIC/, and ftrf -I- hrf + hr^ + nr^ = OmU. By increas- 
ing 71/72, we increase the average number of electrons 
that tunnel through level 1 while level 2 is occupied (the 
quantity Fin in Eq. jZl). Indeed, similar to what we 
expect for the simplified model of Sec. IIIIl (Eqs. ([7]) and 
([5])), the noise is larger for large 71/72 as one can see 
from Fig. ^a). 

In Fig. [71 we plot the left-right asymmetry depen- 
dence of the Fano factor. Similar to the simplified case 
of Sec. mil we find the following asymmetry dependence 
of the Fano factor: for large F^/F^, the Fano factor 
is F = 1, and as we decrease F^/F^, we see enhance- 
ment of the the Fano factor. Also, as we discussed in 
Sec. mil the finite temperature suppresses the Fano fac- 
tor below some F^/F^ ^ 1 toward the value F = 1 
in the limit T'" /T^ — > 0. The most interesting feature 
in Fig. [3 is the unexpected influence that the cotun- 
neling processes have on the Fano factor. At low tem- 
peratures {KbT <C eVbias), the cotunneling processes 
suppress the Fano factor significantly in the asymmet- 
ric coupling regime F^ ^ F''^ (the dashed lines in Fig. [71 
represent calculations without cotunneling processes). 

The physics behind the reduction of the Fano factor in 
the F^/F^ < 1 hmit can be explained by the simplified 
picture of Sec. IIIIl At zero temperature and taking into 
account sequential tunneling only, as we reduce F^/F^ 

the Fano factor is enhanced toward the value F « 7?^ in 

272 

the hmit F^/F^ (see Eqs. ^ and This suits 

the upper (red) dashed line in Fig. [71 Finite temperature 
allows backward tunneling- opposite to the voltage bias 
direction, which reduces the number of current pulses 
through level 1 each time level 2 is occupied. The reason 
is that for F^ <C F^ electrons tunnel many times back 
and forth between the right lead and the system before 
a tunneling event from the left lead to the system takes 
place. Thus, the quantity FiTi in Eq. ([7]) (see also Fig. 
12) is reduced and the Fano factor is suppressed. 
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Figure 7: The Fano factor, F = at £1 = -U (the left 

peak of the current) as a function of the left-right asymme- 
try with the following parameters: eVbtas = 0.1(7, Ei — E2 = 
O.Olf/ and hV^ + ftff + hV^ + hV^ = O.OIC/. We assume that 
the two levels have the same left-right asymmetry, F^/F^. 
When we change p^/F^ we keep the total width of each level 
fixed, (a) The Fano factor for different ratios of the two lev- 
els' widths (7i = h{Vf + rf )/2) with KbT = O.OIC/ and 
eVbias = O.IC/. The dashed lines are the sequential tunnel- 
ing Fano factors and the solid lines are the calculated Fano 
factors including cotunneling processes. In the inside box in 
the upper right corner we zoom in on the small values of 
pLypfl rpj^g Fano factor is larger than one (super-Poissonian 
noise) and maximal in asymmetric coupling, (b) The Fano 
factor as a function of the left-right asymmetry for different 
values of eVbias/r. The ratio between the levels' widths is 



^ 1 



lOr^"^. The dashed lines are the sequential tunneling 



Fano factors and the solid lines are the calculated Fano factors 
including cotunneling processes. The cotunneling processes 
dramatically decrease the Fano factor for large eVbtas/r. 



Cotunneling processes, similar to the finite tempera- 
ture, allow backward tunneling. Virtual tunneling to the 
left is possible as a part of a full two particle cotunneling 
process. The important cotunneling processes which sup- 
press the Fano factor in the asymmetric limit are the in- 
elastic processes than change the system's state between 
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the states |1,0) < — > |0,1). In the asymmetric hmit, the 
total rate of these processes becomes larger than the rate 
of sequential tunneling from the left lead into level 1. As 
a result, the occupation of level 2 changes faster than the 
time delay between consecutive current pulses through 
level 1. The quantity Fiti is reduced and the Fano fac- 
tor is suppressed toward the value _F = 1 in the extremely 
asymmetric limit F^/F^ — >■ 0. 

We want to mention the experimental work that was 
reported in Ref. [2^ , where a super- Poisson noise with a 
strong asymmetry dependence was measured in quantum 
dot that was attached to two leads at strong magnetic 
field. One of the suggested explanations for the enhance- 
ment of the noise in this system was an additional level 
that is weakly coupled to the leads. Indeed, the strong 
dependence of the Fano factor on the asymmetry of the 
dot-leads coupling, which is very similar to the asymme- 
try dependence of the Fano factor we have in Fig. [3 
indicates a possible two level system. We must empha- 
size that although our results fits qualitatively the exper- 
imental measurements^, the rate equations formalism is 
not suitable for quantitative analysis of the experimen- 
tal results since the quantum dot was attached relatively 
strong to the leads, a situation that makes the rate equa- 
tions approach invalid. Nevertheless, qualitatively, we 
believe that the significant asymmetry dependence of the 
Fano factor is a strong evidence for the presence of a sec- 
ond interacting level that was weakly attached to the 
leads in the experimental setup. 

Spinful electrons. In order to avoid unnecessary com- 
plexities, we have assumed that the electrons are spin- 
less. Physically, this situation can be realized by two- 
level quantum dot at strong magnetic field or two single 
level quantum dots at strong magnetic field. The spin 
degree of freedom can be added to the problem in two 
ways: single spinful level with spin dependent couplings 
to the leads or two spinful levels with different couplings 
to the leads. The single spinful level case is very simi- 
lar to the system that we analyze, the only difference is 
the fact that processes with different spins cannot inter- 
fere. Therefore some of the cotunneling rates are slightly 
different {e.g., cotunneling of electron through an empty 
level). Nevertheless, we want to stress that the most im- 
portant cotunneling processes in the more physical situ- 
ation of relatively strong interactions, which are the in- 
elastic processes |1,0) O |0, 1) (or | t) ^ I i))) have 
exactly the same rates as in our model. Moreover, and 
this is the important point, the physics behind the noise 
enhancement in the spinful case and in the spinless case, 
is the same. The noise is enhanced because of the block- 
ing effect which is a result of the Coulomb interactions. 
Therefore, the same qualitative dependence on the asym- 
metry is expected. This is true also in the two spinful 
levels case, and in fact also in multi-level systems with 
interactions. If the blocking effect enhances the noise, we 
expect similar dependence of the Fano factor on /T^- 
with similar suppression due to finite temperature and 
cotunneling. 
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Figure 8: The current and the Fano factor in the small U 
{U < eVbias) region with the following parameters (F = Ff + 

+ T^ + r^): u = 2hr, eFbias = lofir, KbT = nr, Ei - 

E2 = 0.5ftr and 71 = IOOOO72. (a) The current as a function 
of E\ at the symmetric point, F^ = F^. The dashed line 
is the current including sequential tunneling only, and the 
solid line is the current including also cotunneling processes, 
(b) The Fano factor at iJi = as a function of the left- 
right asymmetry. The dashed line is the sequential tunneling 
Fano factor and the solid line is the Fano factor including 
cotunneling processes. The cotunneling processes enhance the 
Fano factor significantly. 



B. Weak interactions 

While U is large, say U > eVbias, the current at the 
left (right) peak, Ei ^ —U {Ei w 0) is changed dra- 
matically when an electron is tunneling into or out of 
level 2, since the strong interaction blocks the tunneling 
through level 1. This is not the case for small U. If 
U is small compared to Vbias the current alternates be- 
tween two relatively close values, and the effect on the 
noise is smaller. Yet, if the ratio between the levels' 
widths, 71/72 is large enough, the enhancement of the 
noise can be important. In Fig. [5] we depict the current 
and the Fano factor for relatively small U {U = eVbias/5) 
with the following parameters (F ee Pf + Pf + Pf + P2 ): 
U = 2hT, eVbias = lOftP, KbT = fiP, Ei - E2 ^ O.bhT 
and 71 = IOOOO72. We see again that the cotunneling 
processes are important as they enhance the Fano factor 
significantly (the cotunneling processes contribute less 
than 20% of the current but almost double the Fano fac- 
tor). Notice also that unlike the Fano factor of the strong 
U case, the Fano factor is maximal when the system is 
symmetrically coupled to the leads, P^ = P^. 
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(a) Sequential tunneling 



(b) Cotunneling 
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Figure 9: The dependence of the tunnehng rates through level 
1 on the occupation of level 2 in the small U (C/ < eVbias) 
region with the following parameters (r = Tf + rf' + r2^+r2 ): 

u = 2hr, Vbias = whv, KbT = hv, Ei = o, E2 = -o.sftr 

and 71 — IOOOO72. Negative cotunneling rates mean that the 
cotunneling processes decrease the current (see for example 
Ref. (3l|l. (a) The sequential tunneling rate from the left lead 
through level 1 to the right lead (((ojj^^oi)"^ + (^oi.ii)""'^)"^ 
when level 2 is occupied and ((tiJio,oo)~^ + (''^00,10)"^)"'^ when 
level 2 is empty), (b) The cotunneling rates from the left lead 
through level 1 to the right lead (wn^n +1^01,01 when level 2 is 
occupied and u^ iq + ojQofio when level 2 is empty), (c) The 
difference between the sequential tunneling rates of the two 
stages (level 2 is empty and full) . (d) The difference between 
the cotunneling rates of the two stages (level 2 is empty and 
full). 



To explain the enhancement of the noise in the small 
U case, we use the intuitive picture that we studied in 
Sec. Illli a sequence of current pulses through level 1 and 
a random telegraph signal describing the occupation of 
level 2. The rate of the pulses depends on the occupation 
of level 2 and we consider two different rates: Fi describes 
the tunneling rate through level 1 while level 2 is full and 
Fi describes the tunneling rate while level 2 is empty. 
The Fano factor is given by ([T^ 



F w 1 



2(AFi) 
TiFa 



(26) 



where the rate F2 ~ {tq + ti) ^ is the tunneling rate 
through level 2, AFi = fi -Fj, Pc = zr^, and we have 
used the fact that f 1 = Fi + AFi 



Ta+Ti ■ 

Fi for [/ < eT4ias- 
We see that although the Fano factor is usually reduced 



to one for ^ 1, in the extreme case of 
the Fano factor can be enhanced. 



We can now understand why the Fano factor is en- 
hanced noticeably due to the cotunneling processes. For 
U < e Vbias, the rate of the sequential tunneling through 
level 1 barely depends on the occupation of level 2 since 
the energies Ei and Ei + U are close to each other com- 
pared to Vbias- The cotunneling processes however, are 
much more sensitive to U. In Fig. [9] we illustrate this 
point by comparing the dependence of the sequential tun- 
neling and the cotunneling rates through level 1 on the 
occupation of level 2, with the same parameters of Fig. 
IHb- In Fig. El^a), we plot the sequential tunneling rates 
from the left lead through level 1 to the right lead when 
level 2 is empty- ((t^I^^^oo)"^ + ('^oo,lo)~^)~^ full- 
((wjj" 01)^"'" + {^01 11)^"'")^"'" (expressions for the rates are 
given in Appendix lAl see Fig. [3] to clarify the notations). 
In Fig. we plot the cotunneling rates from the left lead 
through level 1 when level 2 is empty- lo^ 10 +^^00 oo; ^^^^ 
full- u}]^ II +^^01 01 (expressions for the rates are given in 
Appendix [Aj see Fig. |4] to clarify the notations) . In 
Figs, [nt andinji we plot the difference between the val- 
ues of the tunneling rates through level 1 when level 2 is 
empty and the values of the tunneling rates when level 2 
is full (in Fig. ^ the sequential tunneling rate and|9jl the 
cotunneling rate). We see that although the sequential 
tunneling rate is order of magnitude larger than the co- 
tunneling rate, the cotunneling processes are much more 
sensitive to the occupation of level 2, making the differ- 
ence between the two values of the cotunneling rate on 
the same order of the difference between the two values 
of the sequential tunneling rate. Thus, even though the 
cotunneling processes have a small contribution to the 
current, they contribute the same as the sequential tun- 
neling to AFi and therefore have an important contribu- 
tion to the Fano factor. We want to emphasize that the 
Fano factor is enhanced in the weak interactions regime 
only if level 2 is coupled to the leads extremely weaker 

than level 1, as we require y~ < (t^) where (^t^) 
is a very small number. 

Since (|13p depends quadratically on the multiplication 
of the probabilities of finding level 2 empty and full, 
Po(l — ^"0)1 the enhancement of the Fano factor due to 
AFi is maximal where Pe{l — Pe) is maximal. In the limit 
of small U, the energy of the system, wether it is empty, 
singly occupied or doubly occupies, is more or less the 
same for Ei = 0, therefore Pc(l — ^c) is maximal where 
the system is symmetrically coupled to the leads. 



VI. CONCLUSIONS 

We have analyzed the transport through a two level 
system that is coupled to two leads with one level coupled 
much stronger than the other. We showed that a simple 
intuitive model that describes the transport through this 
system as a sequence of current pulses that is stopped oc- 
casionally captures well most of the transport properties 
of the system. The current pulses in this model corre- 
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spond to tunneling events through the strongly coupled 
level, that are stopped occasionally by tunneling events 
through the weakly coupled level. As expected from this 
simple model, we find in the more rigorous rate equations 
based calculation, a super-Poisson noise, with Fano fac- 
tor larger than one, indicating that the electrons tunnel 
in bunches. 

We showed, using rate equations calculation, a unique 
dependence of the transport on the asymmetry of the 
coupling to the leads, F^/F-^, in the strong interaction 
limit (see Figs [S] and [7]) . In particular, the Fano factor 
is enhanced at asymmetric couplings and increases as we 
decrease F-^/F^. Eventually, the finite temperature sup- 
presses the Fano factor in the limit F^/F^ — > 0. 

Interestingly, at eVbias 3> KbT the cotunneling pro- 
cesses affect dramatically the enhancement of the Fano 
factor. The cotunneling processes that change the sys- 
tem's state between the states |1,0) and |0, 1) (level 1 
occupied and level 2 occupied) suppress the Fano factor 
even for i^sT/eVbias <C 1 and therefore they are essential 
for transport calculations even though most of the cur- 
rent is carried by sequential tunneling processes. Other 
cotunneling processes are less critical in the strong inter- 
action regime, they slightly improve the approximation 
but don't change the transport significantly. 

At relatively weak interactions, we showed that the 
Fano factor may also be enhanced noticeably if the differ- 
ence between the couplings of the two levels is extremely 
large, namely 72/71 ^ 1. In this limit we find again that 
the cotunneling processes are essential for the analysis of 
the noise and the Fano factor. The reason is because the 
telegraphic fluctuations of the sequential processes are 
of the same order as the fluctuations of the cotunneling 
processes. Thus, even though the current is described 
by sequential tunneling quite well, one must include the 
cotunneling processes in the noise calculations. 
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Appendix A: List of tunneling rates 

In this appendix we give the expressions for all the 
tunneling rates that enter the rate equations. 



Sequential tunneling rates 

To the lowest order in Ht the transition rates can be 
calculated using Fermi's golden rule. We use the notation 
for the rate of a tunneling process that changes the 
system's state from 'i' to 'j' by tunneling an electron 
from the left to the right direction (and similarly uj'fj 
for electron that moves from right to left). For example 
Wf)o 10 is the rate of tunneling from the left lead into level 
1 while level 2 is empty. The rates of the sequential 
tunneling processes are 



'00,10 






(Al) 


'00,01 


= F2 Fi?£)(i?2 




(A2) 


'10,00 


= FfFj.,3(M^ 


-El), 


(A3) 


'01,00 


= f^f^,3(m« 


-E2), 


(A4) 


'01,11 


= rf Ff£)(_Ei 




(A5) 


'10,11 


= F2 FF£i(i?2 




(A6) 


'11,01 




-E,-U), 


(A7) 


'11,10 


= r^FFD{f^'' 


-E2-U). 


(A8) 



We have defined 

Ff^.^K^P, Ff^.^lffp, (A9) 

and the Fermi's function Ffd{x) = (1 + e^^)^^ gives the 
probability for the availability of an electron or a hole for 
the tunneling process. The left moving rates, wf^, have 
similar expressions with /i^ ^ jj^ and Ff 2 O Ff 2 • 

Elastic-cotunneling rates 

Each elastic cotunneling process has two possible inter- 
mediate states. For example, electron can tunnel through 
an empty system via level 1 or 2, thus, cotunneling pro- 
cesses of the form |0,0) — > |0,0) have two possible in- 
termediate states: |1,0) and |0, 1). The two possible 
paths interfere and we need to sum the amplitudes of the 
two possibilities rather than their probabilities. We use 
the notation {'-^tl) for elastic-cotunneling processes 
in which the electron tunnels to the right (left) direc- 
tion. The total elastic-cotunneling rates are the sum of 
the rates of all the possible processes, namely integrating 
over all incoming electron's energies. 
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The left moving elastic-cotunneling rates, wf^, have sim- 
ilar expressions with ji^ ^ ji^ . Equations (jA10p - (|A13p 
are formal expressions, and the actual rates, which we use 
as input for the rate equations, cannot be directly calcu- 
lated from these integrals. The reason is the divergence 
of these expressions due to the finite widths of the two 
levels (which we treat as two delta functions in energy). 
We use a regularization scheme^^i^ for the calculation 
of the cotunneling rates. The regularization scheme is 
summarized in Appendix [Bl 



Inelastic-cotunneling rates 

We consider the inelastic-cotunneling processes that 
change the system's state between |1,0) < — > |0, 1) and 

I 



|0,0) i — > |1, 1). We begin with the former; we use the 
notation w,*^ for the rate of processes in which the system 
changes its state from |i) to \ j) in the following way: the 
electron that enters or leaves level 1 tunnels to the right 
direction, while the electron that enters or leaves level 2 
tunnels to the left, and similarly, we define w,^, wf^ . 
For example, if the system's initial state is |1, 0) and the 
electron in level 1 tunnels to the right lead, while another 
electron from the right lead tunnels to level 2, we denote 
the rate of this process by . Again, there are two 

possible intermediate states for the processes 1 1 , 0) < — > 
|0, 1), and we need to sum them properly. The formal 
expression for these rates are 
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(A14) 
(A15) 
(A16) 
(A17) 



To get these integrals we used the energy conservation: if 
the incoming electron (that enters level 2) has the energy 
e, the outgoing electron (that leaves level 1) must have 
the energy e + Ei — £2- To get the rates: u)^ ^q, uj^ ^g, 
'^01 10 ^^^"^ "^01 10' may write integrals like (jA14p - (|A17|) 

L R L R 

and exchange tj^ ' 0^2' ^md Ei -(-> £2 • Once again, the 
formal expressions (jA14|) - (|A17p need to be regularized in 
order to extract the input terms for the rate equations 
calculation (see Appendix [BJ. 



The last rates that we discuss are the inelastic 
|0,0) i — > processes' rates. These processes are 

somewhat more complex than the other cotunneling pro- 
cesses as they have four possible intermediate states. We 
use the notation n (w^^^) and qq ("^uoo) 
processes in which two electrons enter or leave the two 
levels by tunneling to the right (left). We use the no- 
tation Wqq II, QQ for processes in which the two elec- 
trons enter or leave the two levels by tunneling one to 
the right and the other to the left. The rates for the 
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|1, 1) inelastic cotunneling processes are 
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The factor 1/2 that appears in Eqs. (|AT8|) . ([AT9|) . 
(jA2ip and (|A22p is due to the double counting of pro- 
cesses: by integrating over e we sum both e = e' and 
e = — e' + + i?2 + however these two processes 
are identical since in both cases when the two levels are 
empty the two electrons are in the same lead with ener- 
gies e' and —e' + Ei+ E2 + U. Hence, as we double count 
each process we insert a factor of 1/2. Notice also that 
for t^TO.ii' '^ooai' '^n.oo ^^'^ '^n.oo' by including the point 
e — ^^+^^+^ in the integral, we include an impossible 
process, as the two electrons in the lead have the same en- 
ergy Nevertheless, this point contribute zero 
to the integrals, and therefore we have no problems with 
the formal expressions (|Jl8|, ((JT9|l . (|X2T|) and (|J22|) . 
The integrals in (Pl8)) - (|M3| need a regularization in or- 
der to extract the rates that we use in the rate equations, 
the regularization scheme appears in Appendix |B] 



where the amplitudes A, B are multiplications of two tun- 
neling coefficients, and one might need to use the relation 
Ffd{~^) = ^ — Ffd{() hi order to bring the expression of 
a specific cotunneling rate to this form. The integral (jBl[) 
diverges due to the finite widths of the energy levels^^iS. 
We first add by hand a width to the levels 



I{A,B,Ea,Eb,fIl,ll2) = 



h 



(B2) 
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Appendix B: APPENDIX: REGULARIZATION 
SCHEME 



In this appendix we summarize the regularization pro- 
cedure for the cotunneling rates^^i^. AH the cotun- 
neling rates that appear in Appendix |^ excluding the 
|0,0) i — > |1, 1) rates, can be brought to the form 



I{A,B,Ea,Et,^ll,^i2) = X 



/ 



h 



deFpoie - Mi)(l - Fpoie - ^2)) 



(Bl) 



A B 

+ 



e-Ea e-Et 



Next, we solve the integral (jB2p and write the solution 
as a power series in 7. We extract the cotunneling rate by 
subtracting the I/7 term and taking the limit 7 — > 0. We 
should emphasize that in general each level has its own 
width and the sign ±17 depends on the process; we asso- 
ciate difl^erent signs for incoming and outgoing electrons. 
If one adds the widths properly, the divergent term (I/7) 
has a physical meaning, and one can read the sequential 
tunneling rates from it. Nevertheless, these details are 
not important for the regularization procedure, and the 
finite values of the cotunneling rates are independent of 
the details of the regulator 7. We can write (jB2p as 



I{A,B,Ea,Eb,^ii,lJL2) 



2^z/2 



(n(A, Ea, III, 112) + Il{B, Eb, ni,H2) + /2(A, B, Ea, Eb, A*i, A/2)) 
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where 
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The solutions of II and 12 can be written using digamma functions with complex variables (iplz)). II contains a 
divergent part 
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where /3 = 1/{KbT) and ^s(M2 - = (e'^^^^^^i^ - 1)"^ After the subtraction of the I/7 term and taking the 
limit 7 — > 0: 
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There is no divergence in 12, therefore we simply solve it and take the limit 7 — > 
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We can use the solution of (jBip to solve the rates of to write Eqs. (|A18P - (jA23|) as sums of six terms of the 
the processes |0, 0) < — >■ |1,1). Notice that the formal form of (jBip . 
expressions for these rates (Eqs. (jA18P - (jA23P ') contain 
four terms inside the absolute value. We may use the 
trivial identity 
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